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We study synchronization in the context of network traffic on a 2 — d communication network with 
local clustering and geographic separations. The network consists of nodes and randomly distributed 
hubs where the top five hubs ranked according to their coefficient of betweenness centrality (CBC) are 
connected by random assortative and gradient mechanisms. For multiple message traffic, messages 
can trap at the high CBC hubs, and congestion can build up on the network with long queues at the 
congested hubs. The queue lengths are seen to synchronize in the congested phase. Both complete 
and phase synchronization is seen, between pairs of hubs. In the decongested phase, the pairs start 
clearing, and synchronization is lost. A cascading master-slave relation is seen between the hubs, 
with the slower hubs (which are slow to decongest) driving the faster ones. These are usually the hubs 
of high CBC. Similar results are seen for traffic of constant density. Total synchronization between 
the hubs of high CBC is also seen in the congested regime. Similar behavior is seen for traffic on 
a network constructed using the Waxman random topology generator. We also demonstrate the 
existence of phase synchronization in real Internet traffic data. 

PACS numbers: 89.75.Hc 



6 
o 

O 
• I— I 



> 

p 
o 

00 

o 



X 



I. INTRODUCTION 



The phenomenon of synchronization has been studied 
in contexts ranging from the synchronization of clocks 
and the flashing of fire-flies [l| to synchronization in os- 
cillator networks 0] and in complex networks 0]. Syn- 
chronized states have been seen in the context of traf- 
fic flows as well Q, and investigations of traffic fiow on 
substrates of various geometries have been the focus of 
recent research interest [1, 0, S B ■ The synchronization 
of processes at the nodes, or hubs, of complex networks 
can have serious consequences for the performance of the 
network Q. In the case of communication networks, the 
performance of the networks is assessed in terms of their 
efficiency at packet delivery. Such networks can show 
a congestion-decongestion transition [13. We note that 
an intimate connection between congestion and synchro- 
nization effects has been seen in the case of real networks 

The aim of this paper is to study the interplay of con- 
gestion and synchronization effects on each other, and 
examine their effect on the efficiency of the network for 
packet delivery in the context of two model networks 
based on two dimensional grids. The first network con- 
sists of nodes and hubs, with the hubs being connected 
by random assortative or gradient connections [l3| . In the 
case of the second network, in addition to nearest neigh- 
bour connections between nodes, the nodes are connected 
probabilistically to other nodes, with the probability of 
a connection between nodes being dependent on the Eu- 
clidean distance between them[lj|. Such networks are 
called Waxman networks and are popular models of in- 
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ternct topology [l5j. Synchronisation effects are observed 
in the congested phase of both these model networks. In 
addition to these two networks, we also discuss synchro- 
nisation effects seen in actual internet data. 

We first study synchronization behavior in a two di- 
mensional communication network of nodes and hubs. 
Such networks have been considered earlier in the con- 
text of search algorithms Il6| and of network traffic with 
routers and hosts [13, [ll,ll9[. Despite the regular 2 — d 
geometry such models have shown log-normal distribu- 
tion in latency times as seen in Internet dynamics [20| . 
The lattice consists of two types of nodes, the regular 
or ordinary nodes, which arc connected to each of their 
nearest neighbors, and the hubs, which are connected 
to all the nodes in a given area of influence, and are 
randomly distributed in the lattice. Thus, the network 
represents a model with local clustering and geographi- 
cal separations [2l|, [l^. Congestion effects are seen on 
this network when a large number of messages travel be- 
tween multiple sources and targets due to various factors 
like capacity, band- width and network topology [23| . De- 
congestion strategies, which involve the manipulation of 
factors like capacity and connectivity have been set up for 
these networks. Effective connectivity strategies have fo- 
cused on setting up random assortative [2^, or gradient 
connections |25j between hubs of high betweenness cen- 
trality. 

We introduce the ideas of phase synchronization and 
complete synchronization in the context of the queue 
lengths at the hubs. The queue at a given hub is de- 
fined to be the number of messages which have the hub 
as a temporary target. During multiple message transfer, 
when many messages run simultaneously on the lattice, 
the network tends to congest when the number of mes- 
sages exceed a certain critical number, and the queue 
lengths tend to build up at hubs which see heavy traffic. 
The hubs which see heavy traffic are ranked by the co- 
efficient of betweenness centrality {CBC), which is the 
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fraction of messages which pass through a given hub. 
We focus on the top five hubs ranked by CBC. Phase 
synchronization is seen between pairs of hubs of compa- 
rable betweenness centrahty. The hub which is slowest 
to decongest (generally the hub of highest CBC) drives 
the slower hubs with a cascading master-slave effect in 
the hub hierarchy. When the network starts dccongcst- 
ing, the queue lengths decrease, and synchronization is 
lost. These results are reflected in the global synchroniza- 
tion parameter. When decongestion strategies which set 
up random assortative, or gradient, connections between 
hubs are implemented, complete synchronization is seen 
between some pairs of these hubs in the congested phase, 
and phase synchronization is seen between others. We 
demonstrate our results in the context of the gradient 
decongestion strategy, but the results remain unaltered 
for decongestion strategics based on random assortative 
connections. Similar results are seen for constant density 
traffic where a fixed number of messages are fed on the 
system at regular intervals. Total synchronization is also 
seen in the queue lengths of the hubs of high CBC. 

All the results obtained for the first model are ob- 
served for message transport on Waxman topology net- 
work, where again synchronisation of hubs of high CBC 
is observed in the congested state. We demonstrate these 
results. Finally we study internet traffic data and demon- 
strate that phase synchronisation is seen in this data as 
well. Intermittent phase synchronization is also seen in 
this data. 



II. A COMMUNICATION NETWORK WITH 
LOCAL CLUSTERING AND GEOGRAPHIC 
SEPARATION 

We first study traffic congestion for a model network 
with local clustering developed in Ref. [l^ . This network 
consists of a two-dimensional lattice with ordinary nodes 
and hubs (See Fig. [1]). Each ordinary node is connected 
to its nearest- neighbors, whereas the hubs are connected 
to all nodes within a given area of infiuence defined as 
a square of side 2k centered around the hub[T3|. The 
hubs are randomly distributed on the lattice such that 
no two hubs are separated by less than a minimum dis- 
tance, dmin- Constituent nodes in the overlap areas of 
hubs acquire connections to all the hubs whose influence 
areas overlap. The source S{is,js) and target T(it,jt) 
are chosen from the lattice and separated by a flxed dis- 
tance Dst which is defined by the Manhattan distance 
Dst = \is — it\ + \js — jt\ . It is useful to identify and 
rank hubs which see the maximum traffic. This is done by 
defining the co-efficient of betweenness centrality (CBC) 
where the CBC of a given hub k is defined as CBC = 
i.e. the ratio of the number of messages that go through 
a hub k to the total number of messages running on the 
lattice. These are listed in Table HI 

Efficient decongestion strategies have been set up by 
connecting hubs of high CBC amongst themselves, or to 



randomly chosen other hubs via assortative connections 
[23 |. Gradient mechanisms [l^] can also be used to de- 
congest traffic [2l[23(See Fig. [^b)). 

In all the simulations here, we consider a lattice of size 
100 X 100 with 4% hub density and D^t = 142, d„„ = 1. 
The critical message density which congests this lattice 
is Nc ~ 1530. The studies carried out here correspond 
to the congested phase, where 2000 or 4000 messages run 
on the lattice. We first consider the baseline lattice as 
in Fig. [ija) where there are no short-cuts between the 
hubs. The message holding capacity of ordinary nodes 
and hubs is unity for the baseline lattice. 

A given number N of source and target pairs sepa- 
rated by a fixed distance Dst are randomly selected on 
the lattice. Here, all source nodes start sending messages 
to the selected recipient nodes simultaneously, however, 
each node can act as a source for only one message during 
a given run. The routing takes place by a distance based 
algorithm in which each node holding a message directed 
towards a target tries to identify the hub nearest to it- 
self, and in the direction of the target as the temporary 
target, and tries to send the message to the temporary 
target through the connections available to it. During 
peak traffic, when many messages run, some of the hubs, 
which are located such that many paths pass through 
them, have to handle more messages than they are capa- 
ble of holding simultaneously. Messages tend to jam in 
the vicinity of such hubs (usually the hubs of high CBC) 
leading to formation of transport traps which leads to 
congestion in the network. Other factors like the oppos- 
ing movement of messages from sources and targets situ- 
ated on different sides of the lattice, as well as edge effects 
ultimately result in the formation of transport traps. We 
have studied trapping configurations for the same 2 — d 
network in Ref. [2^ . Fig. UKa) shows a situation in which 
messages are trapped in the vicinity of high CBC hubs. 
Fig. Wijo) shows the number of messages running on the 
lattice as a function of time. It is clear that the messages 
are trapped for the baseline case. We study the network 
for situations which show this congested phase. 



III. QUEUE LENGTHS AND 
SYNCHRONIZATION 

As mentioned in the introduction, the queue at a given 
hub is defined to be the number of messages which have 
the hub as a temporary target. As traffic increases in 
the network, hubs which see heavy traffic start getting 
choked due to capacity limitations, and are unable to 
transfer messages aimed towards them to the next tem- 
porary target. Thus queue lengths start to build up at 
these hubs. If these hubs arc not decongcstcd quickly, so 
that the queue lengths start falling, the congestion starts 
spilling over to other hubs. If the number of messages 
increases beyond a certain critical number, messages get 
trapped irretrievably and the entire lattice congests. A 
plot of the queue lengths as a function of time can be seen 
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TABLE I: This table shows the CBC values and ranking of the top five hubs. A total number of 2000 messages are traveling 
simultaneously on a 100 x 100 lattice with 4% hub density and Dst = 142 and run time set at 5000. 

Hub label CBC value Rank 

^ 0827 1 

y 0.734 2 

z 0.726 3 

u 0.707 4 

V 0.705 5 




FIG. 1; (Color online) (a) A regular two dimensional lattice. The node X\s an ordinary node with nearest neighbor connections. 
Each hub has a square influence region (as shown for the hub Y). A typical path from the source S to the target T is shown with 
labeled sites. The path 5- 1 -2-3-A-4-5-6- 7-B-8-9- lO-T passes through the hubs A and B. (b) The high CBC hubs 
(M-Q) are connected by the CBCa mechanism. A message is routed along the shortest path S\ — a— M— N— O— 6 — c — d — T\. 
We enhance the capacities of high CBC hubs (A-G) proportional to their CBC values by a factor of 10. Connections between 
these hubs are made by the gradient mechanism. After the implementation of the gradient mechanism, the distance between 
G and F is covered in one step as shown by the link g and a message is routed along the path S — 1 — 2 — 3— G— g— F— 4 — 5 — 6 — T. 



in Fig. [3]Ja). Here, the queues at the 1st, '^nd and 5th hubs 
ranked by the CBC are plotted for the base-hne network 
with no decongcstion strategies implemented. Thus the 
network congests very easily. Since the queue length is 
defined as the number of messages with the given hub 
as the temporary target, the queue starts dropping as 
soon as the hub starts clearing messages and reaches a 
minimum. Meanwhile, other hubs which were temporary 
targets have cleared their messages, and some new mes- 
sages pick up the hub of interest as their temporary tar- 
get. The queue thus starts building up here, and reaches 
a maximum. After this the messages start clearing, and 
the queues drop sharply. However, since the number of 
messages is sufficiently large for the network to congest, 
some messages get trapped in the vicinity of the hub, and 
the queues saturate to a constant value. Similar phenom- 
ena can be seen at the hubs of lower CBC (see Fig. [Sja)). 
Here again three distinct scales can be seen with values 
of the same order as those for the highest ranked hub. 



An important difference can be seen in the queues of the 
fifth ranked hub (Fig. ^a)) as well as the fourth ranked 
hub (not shown). Since these hubs have lower CBC val- 
ues, and thus fewer messages take them as the temporary 
targets, the queues at these hubs clear completely. Thus, 
the saturation value at these hubs, is zero. It should 
also be noted that the time at which the last two hubs 
clear completely, i.e. the queue length drops to zero, is 
substantially earlier than the saturation time of the top 
two hubs. The above results are observed for a typical 
configuration and is valid for different configurations as 
well. 



A. Synchronization 

We now study the synchronization between the queues 
at different hubs. We see phase synchronization between 
queues at pairs of high CBC hubs for the baseline, and 
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(a) 1 (b) t 

FIG. 2: (Color online) (a) Figure shows messages trapped in the vicinity of top five hubs ranked by their CBC values. Numbers 
indicate the ranks of the top five hubs. A total number of 133 messages are trapped in the lattice when 2000 messages are 
traveling simultaneously on a 100 x 100 lattice with 4% hub density for a run time of lODst- (b) Messages are trapped in the 
lattice when the hubs are not connected (baseline mechanism). All messages are delivered to their respective targets once the 
top five hubs ranked by their CBC values are connected by the gradient mechanism. 



complete synchronization between some pairs once de- 
congestion strategics arc implemented. The usual defi- 
nitions of complete synchronization and phase synchro- 
nization in the literature are as follows. 

Complete synchronization (CS) in coupled identical 
systems appears as the equality of the state variables 
while evolving in time. Other names were given in the lit- 
erature, such as conventional synchronization or identical 
synchronization^^ . It has been observed that for chaotic 
oscillators starting from uncoupled non-synchronized os- 
cillatory systems, with the increase of coupling strength, 
a weak degree of synchronization, the phase synchro- 
nization(PS) where the phases become locked is seen 
[sol . [3l| . Classically, the phase synchronization of cou- 
pled periodic oscillators is defined as the locking of phases 
01,2 with a ratio n : m (n and m are integers), i.e. 
1 7101 ~ '71021 < Const. 

These two concepts of synchronization are applied to 
the queue lengths of the top five hubs. The plot of qi{t) 
as a function of average queue length < q{t) > shows 
a loop in the congested phase, similar to that observed 
in coupled chaotic oscillators [331. We define a phase as 
in [s^l ^i{t) = t^'iT-^^ { ^q(t)> )' "^here qi{t) is the queue 
length of ith hub at time t, and < q{t) >= qj (t) 

i 

where the average is calculated over the top five hubs 
{Nfi = 5). The queue lengths are phase synchronized if 

\<P,{t) - <i>j{t)\ < Const (1) 

where ^i{t) and $j(t) are the phase at time t of the ith 
and jth hub respectively. 



Two queue lengths qi{t) and qj{t) arc said to be com- 
pletely synchronized if 

q^it)^q,it) (2) 

Fig. El^b) shows that the queue lengths of the first 
and fifth ranked hubs are not completely synchronized. 
Fig. [3]Jc) shows the phase difference between the top 
pair of hubs as a function of time for the base-line case. 
It is clear that the two hubs are phase synchronized in 
the regimes where the queues congest. There are three 
distinct time scales in the problem. The two hubs are 
phase synchronized up to the first time scale ti , where the 
queues cross each other first, they lose synchronization 
after this. The point at which the phase difference is 
maximum is t2. This is the point at which the first hub 
saturates, but the second hub is still capable of clearing 
its queue. At ^3 both the hubs get trapped and the phases 
lock again. 

Fig. Eljd) shows a similar plot for the hubs of the 
two remaining ranks. It's clear that the hubs phase syn- 
chronize. The synchronization behavior of the remaining 
hubs for a typical configuration is listed in Table [TTl It 
is clear that the hubs synchronize pair wise, and that 
the slower hubs drive the hubs which clear faster. Since 
the queues at the fourth and fifth hub clear faster than 
the first hub saturates, there is no peak in the A^jisj 
plot for the (1, 5) pair and hence no scale ^2- The phase 
synchronization between hubs three and four shows sim- 
ilar behavior. This is valid for different configurations as 
well. 

It is also interesting to compare the synchronization 
effects between these hubs of high CBC, and randomly 
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FIG. 3: (Color online)We consider the baseline mechanism, (a) Queue lengths for the hubs ranked 1, 2 and 5 ranked by 
CBC. (b) Complete synchronization is not observed in the queue lengths. The error in the fit to j/ = a; is 1.940. (c) Phase 
synchronisation of the queue lengths of the two top most hubs, (d) Phase synchronization of the queue lengths of the 1st and 
5th hub ranked by CBC. A total number of 2000 messages are traveling simultaneously on a 100 x 100 lattice with 4% hub density. 



selected hubs on the lattice. Fig. U] shows the phase dif- 
ference between the hub of highest CBC (hub 'x', ranked 
1) and a randomly chosen hub (with CBC value 0.56 ). 
It is clear that there is no synchronization between these 
two hubs. However, this randomly chosen hub shows 
excellent phase synchronization with another randomly 
chosen hub (with comparable CBC value 0.67). Similar 
results are seen for larger number of messages. 



B. Decongestion strategies and the role of 
connections 



As discussed earlier, the addition of extra connections 
between the hubs of high betweenness centrality can ease 
congestion. This leads to two effects. The time scales of 
the problem, the rate at which the queues build up and 
clear, and the way in which correlations occur between 
different hubs are altered due to the addition of extra 
connections. We see the effects of this in the synchro- 
nization between the queues at the hubs. We illustrate 
the effects seen for gradient connections between the hubs 
(Fig. mb)). To set up the gradient mechanism, wc en- 
hance the capacities of top five hubs ranked by their CBC 
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TABLE II: The table shows the pair of queue lengths of the top five hubs which are phase synchronized for a typical configuration 
for the baseline mechanism. All the pairs of top five hubs are phase synchronized within a constant C — 0.05. The three time 
scales ti, t2 and for every synchronized pair are shown in the table below. Similar results are observed for different 
configurations as well. We consider a 100 x 100 lattice with 4% hub density and D^t ~ 142 . A total number of 2000 messages 
are traveling simultaneously in the lattice. The run time is set at 5000. 
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FIG. 4: (Color online) No phase synchronization is observed between top most hub and a randomly chosen hub. The two 
randomly chosen hubs are phase synchronized to large extent. All parameters as in Fig. [S] A total number of 2000 messages 
are traveling simultaneously on a 100 x 100 lattice with 4% hub density. 



values, proportional to their CBC values by a factor of 
10. A gradient flow is assigned from each hub to all 
the hubs with the maximum capacity (Cmax)- Thus, the 
hubs with lower capacities are connected to the hubs with 
highest capacity Cmax by the gradient mechanism. Hence 
the hub with highest CBC value is maximally connected. 
Fig. ^b) and Fig. [5] shows that connecting the top five 
hubs by the gradient mechanism relieves the system of 
congestion rapidly when 2000 messages are traveling in 
the lattice for 4% hub density and run time of lODgt- 

The most striking observation is that now complete 
synchronization is seen between at least one pair of hubs, 
and phase synchronization is seen between the remain- 
ing pairs. In Fig. ^a.) we plot a pair of queue length 
Qi vs Qj. If these two quantities lie along the y = x line 
with a standard deviation less than one, we call them 



completely synchronized. It is clear from the Fig. [U^a) 
and the value of the standard deviation that the queue 
lengths of the 2nd and the 3rd ranked hubs are completely 
synchronized. These two hubs are of comparable CBC 
values (See Table H]) and are indirectly connected via the 
top most ranked hub, to which each of the lower ranked 
hubs is connected via a gradient. If the standard de- 
viation is greater than one, the queue lengths are not 
completely synchronized. In Fig. [UJb) we observe phase 
synchronization when the top five hubs are connected by 
the gradient mechanism. Phase synchronization is ob- 
served when the queues congest. As soon as the queues 
decongest they are no longer phase synchronized. This 
observation is true for the complete synchronization as 
well. In the gradient scheme we see a star-like geometry 
where the central hub is connected to the hubs of low 
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FIG. 5: (Color online) Plot of queue lengths of the 1st, 2„£j and Ath hubs for the gradient mechanism. A total number of 2000 
messages are traveling simultaneously on a 100 x 100 lattice with 4% hub density. 




(a) ^3 (b) t 

FIG. 6: (Color online) Figure shows (a) complete synchronization in queue lengths for the 2„d and 3rd hub and (b) phase 
synchronization in queue lengths for pair(s) of top five hubs when connected by the gradient mechanism. A total number 
of 2000 messages are traveling simultaneously on a 100 x 100 lattice with 4% hub density. The error in the fit y = a:: in (a) is 0.992. 



capacity. This central hub gets congested leading to the 
congestion of the rest of the hubs. Once this hub gets de- 
congested the rest of the hubs of high CBC get cleared. 
Thus, the central hub, which is the hub of highest CBC, 
drives the rest. 



C. The finite time Lyapunov exponent 

The queue lengths increase in the congested phase and 
the difference between two queue lengths (i.e. queue 
lengths at distinct hubs) is small in this phase, as com- 
pared to the decongested phase, where the difference be- 
tween queue lengths is large. This is analogous to the 
behaviour of trajectories in the chaotic regime where 
the separation between two co-evolving trajectories with 
neighbouring initial conditions increases rapidly, as com- 
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pared to the separation in the periodic regime where it 
rapidly decreases. Hence the stabiHty of the completely 
synchronized state seen in the gradient case can studied 
by calculating the finite time Lyapunov exponent of the 
separation of queue lengths for the top five pairs of hubs. 
The finite time Lyapunov exponent is given by 



AW = iln(M) 
^' t ^5{0)' 



(3) 



where 5{t) — \qi{t) ~Qj{t) \ and (5(0) is the initial differ- 
ence in queue lengths [40|. If \{t) < then queue lengths 
are completely synchronized (Fig. [7Ja)) and if \{t) > 
then queue lengths are not completely synchronized (Fig. 
[71(b)). The time is calculated from the time {t = 15) at 
which the queue starts building up in the lattice. It is 
clear from the Fig. El^a) that complete synchronization 
exists till tc = 720. This is the time at which queues are 
cleared. In Fig. [7Kb) complete synchronization exists till 
tcs = 150; when queues are building up in the lattice. 
No complete synchronization is observed after this, but 
queues are cleared at = 740 . 



D. Global synchronization 



and CBCc (two way)) and the top 5 hubs and ran- 
domly chosen other hubs (Ci3Cb(one way) and CBCd 
(two way)) have also proved to be effective. The phe- 
nomena of complete synchronization and phase synchro- 
nization can be seen for these schemes as well. (See Ta- 
ble imp . Apart from the gradient mechanism complete 
synchronization is seen for the CBCc mechanism as well, 
where the Ath and ^th ranked hubs (ranked by CBC) 
are completely synchronized. Unlike the gradient mech- 
anism, the Ath and bth ranked hubs have a direct two 
way connection for this realization of the CBCc mech- 
anism. Both these hubs have comparable CBC values 
(See Table[I]) and therefore, we see that the queue lengths 
arc completely synchronized. The error to the fit to the 
y = X line is 0.926. The FTLE of the queue lengths 
of these hubs is less than zero indicating complete syn- 
chronization. No complete synchronization is observed 
for the other assortative mechanisms. Its clear from Ta- 
ble [TIT] that the top most hub (labeled x) drives the rest 
of the top five hubs. Global synchronization emerges for 
these cases as well. Thus it is seen that synchronization 
in queue lengths is a robust phenomena. Irrespective of 
the nature of connections between high CBC hubs, syn- 
chronization in queue lengths of highly congested hubs 
exists during the congested phase. 



It is useful to define an over-all characterizer of emerg- 
ing collective behavior. The usual characterizer of global 
synchronization is the order parameter (ssj defined by 

1 ^ 

r exp iip = exp i^j (4) 

iV/i=5, where we consider the top five hubs. Here ip rep- 
resents the average phase of the system, and the $j-s 
are the phases defined in Eq. 1. Here the parameter 
< r < 1 represents the order parameter of the system 
with the value r = 1 being the indicator of total synchro- 
nization. 

We plot the order parameter r and average phase ip 
as a function of time for the baseline mechanism in Fig. 
[SJa) and the gradient connections in Fig. [5]Jb). It is 
clear that the order parameter goes to one up to the 
time tsynch — 300 indicating that the queues at all the 
top 5 hubs synchronize up to this point. As discussed 
earlier, this point is also the time at which the network 
congests. Thus the intimate connection between conges- 
tion and synchronization is clearly demonstrated by the 
order parameter [47|. It is to be noted that the order pa- 
rameter r and average phase ip are calculated for hubs of 
comparable CBC values (in this case the top five hubs). 

E. Other decongestion schemes 

Decongestion schemes based on random assortative 
connections between the top five hubs (CSCq (one way) 



IV. A NETWORK WITH RANDOM WAXMAN 
TOPOLOGY 

The random network topology generator introduced by 
Waxman [31 is a geographic model for the growth of a 
computer network. In this model the nodes of the net- 
work are uniformly distributed in the plane and edges are 
added according to probabilities that depend on the dis- 
tances between the nodes. Such networks are useful for 
Internet modeling due to the distance dependence in link 
formation which is characteristic of real world networks 
[ist and have been widely used to model the topology of 
intra-domain networks [s^, H^, [H, [IJ . We study queue 
length synchronization on this network, and compare it 
with the synchronization seen for the geographically clus- 
tered network of Section II. We consider the case where 
the Waxman graphs are generated on a rectangular co- 
ordinate grid of side L with a probability P(a, b) of an 
edge from node a to node b given by 

P(a,6)=/3exp(--^) (5) 

where the parameters < a, (3 < 1, d is the Euclidean 
distance from a to 6 and M = \f2 x L is the maximum 
distance between any two nodes [H, [H, [3§| . Larger val- 
ues of /? results in graphs with larger link densities and 
smaller values of a increase the density of short links as 
compared to the longer ones. 

Here we select a 100 x 100 lattice. A topology similar 
to Waxman graphs is generated by selecting randomly a 
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FIG. 7: (Color online) Plot of the Lyapunov exponent for (a) complete synchronization in queue lengths of 2nd and 3rd hub 
and (b) No complete synchronization in queue lengths of 1st and 5th hubs when connected by the gradient mechanism, tc is 
the time at which all the messages are delivered to the respective targets. In (a) tc = 720 and (b) tc = 740. A total number of 
2000 messages are traveling simultaneously on a 100 x 100 lattice with 4% hub density. 
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FIG. 8: (Color online) Plot of order parameter r and average phase i/; as a function of time for (a) the baseline mechanism 
and (b) the gradient mechanism. A total number of 2000 messages are traveling simultaneously on a 100 x 100 lattice with 
4% hub density. 



pre-determined number of nodes in the lattice. The 
nodes are then connected by the Waxman algorithm, re- 
sulting in a topology which is similar to Waxman graphs 
(Fig. [5]). Additionally, each node has a connection to 
its nearest neighbors. We study message transfer by the 
same routing algorithm as used in Section II. We evalu- 
ate the coefBcient of betweenness centrality of the nodes 
and select the five top most nodes ranked by their CBC 



values. We compare the synchronization in queue lengths 
of these nodes for different values of (3 and a for simul- 
taneous message transfer. 

The phenomenon of phase synchronization in queue 
lengths is again studied for simultaneous message transfer 
where N = 2000 messages fiow simultaneously on the 
lattice with ~ 100 points chosen randomly in the 
lattice and connected by Waxman algorithm. The source 
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TABLE III: The table shows the pair of queue lengths of the top five hubs which are completely synchronized (with errors to 
the fit to the y — x line in []-s) and phase synchronized for a typical configuration. Similar results are obtained for different 
configurations as well. We consider a 100 x 100 lattice with 4% hub density and Dst ~ 142 . A total number of 2000 messages 
are traveling simultaneously in the lattice. The run time is set at 4:Dst . 



Mechanism 



Complete Synchronization 



Phase Synchronization 



CBCa 

CBCb 
CBCc 
CBCd 
Gradient 



(u,v) [0.926] 
(y,z) [0.992] 



(x,y),(x,z),(x,u),(x,v),(y,z),(y,u),(y,v),(z,u),(z,v),(u,v) 
(x,y) , (x,z) , (x,u) , (x,v) , (y,z) , (y,u) , (y,v) , (z,u) , (z,v) , (u,v) 

(x,y),(x,z),(x,u),(x,v),(y,z),(y,u),(y,v),(z,u),(z,v) 
(x,y),(x,z),(x,u),(x,v),(y,z),(y,u),(y,v),(z,u),(z,v),(u,v) 

(x,y),(x,z),(x,u),(x,v),(y,u),(y,v),(z,u),(z,v),(u,v) 




FIG. 9: (Color online) Figure shows a Waxman topology network generated by connecting 55 points by the Waxman algorithm 
for a — 0.05 and /3 — 0.1, on a 10 x 10 lattice. The number of links increases as the values of a and (3 are increased. 



target separation is Dst = 142 as before, and is again 
the Manhattan distance between source and target. If 
a = 0.05 and /3 = 0.05, the number of Unks between 
the randomly distributed nodes are very few as in the 
geographically clustered network. In such a situation, 
messages are cleared slowly and we observe strong phase 
synchronization (Fig. I10( a)(i)). An increment in the 
values of a and /? increases the density of links. Messages 
are cleared faster due to the presence of a large number of 
short cuts which leads to larger fluctuations in phase, and 
weaker phase synchronization is seen (Fig. [IW a^di)). For 
both the situations messages get trapped and after some 
time and the phase gets locked. We see the saturation 
in the plot of A<I>. Global synchronization is also seen in 
this system (Fig. ^hj). 



V. CONSTANT DENSITY TRAFFIC 

In the previous sections we discussed synchronization 
in queue lengths for simultaneous message transfer where 



messages arc deposited simultaneously on the lattice 
and no further messages are fed on to the system. In this 
section we study synchronization in queue lengths for the 
constant density traffic. For the model of section II we 
consider 100 messages fed at every 120 time steps with 
100 hubs and Dst = 142 for a total run time of 12000. 
Again, two phases, viz. the decongested phase and the 
congested phase are seen. In the decongested phase all 
messages are delivered to their respective targets, despite 
the fact that new messages are coming in at regular in- 
tervals. The queue lengths are not phase synchronized 
during this phase. In the congested phase, messages tend 
to get trapped in the vicinity of the hubs of high CBC, 
due to the reasons discussed in Section II. As more mes- 
sages come in, the number of undelivered messages in- 
crease and the queue lengths start increasing until total 
trapping occurs in the system. During this phase, the 
creation of messages is stopped and the system attains 
maximal congestion. The queue lengths show phase syn- 
chronization during this phase (Fig. [TlT a)). Initially the 
fluctuations in A$ are large. After a time t ~ 2000 the 
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(a) t (b) t 

FIG. 10: (Color online) (a) Plot of phase synchronization in queue length of the top most congested nodes ranked by their CBC 
values for (i) a = 0.05 and (3 — 0.05; (ii) a — 0.4 and /? — 0.4 for simultaneous message transfer. A^™ = 100 points are chosen 
randomly in the lattice and connected by Waxman algorithm. We drop A*" = 2000 messages simultaneously on the lattice, (b) 
The plot of global synchronization parameter r and 'I' as a function of time t for a — 0.05 and /3 — 0.05 ; a = 0.4 and (3 — 0.4 
for simultaneous message transfer. 



queue lengths start increasing and the fluctuations are 
reduced indicating stronger phase synchronization. As 
soon as maximal congestion takes place (~ 7800), the 
phase difference attains a constant value. Global syn- 
chronization is also seen in this system as can be seen 
from Fig. [TlTb) . Note that the scales on which the phase 
difference and the global synchronization parameter fluc- 
tuate is very small indicating a much stronger version of 
synchronization than in the earlier case. 

The results are compared with constant density traffic 
for the Waxman topology network as discussed earlier. 
We consider 100 messages fed continuously at every 120 
time steps for a total run time of 50000 with = 500 
points and Dst = 142. If a = 0.05 and f3 = 0.05 messages 
get trapped in the system very fast {tc = 8000). Phase 
synchronization in queue lengths is observed in such cases 
(Fig. [TlTc)fi)). If the values of a and (3 increase the num- 
ber of links increase. Phase synchronization in queue 
lengths take place at much higher time for (i) a = 0.05 
and P = 0.4 and (ii) a = 0.4 and /3 = 0.05. If a = 0.4 
and P = 0.4 the density of links is very large and all the 
top flve nodes have approximately equal queue lengths. 
Hence we observe a stronger phase synchronization where 
the fluctuation of A$ is well below the predetermined 
constant C = 0.05 (Fig. [TlT c)(ii)). Global synchroniza- 
tion is also observed in this system (Fig. [TTT d)). 

Thus, the model networks studied here show phase syn- 
chronisation as well as global synchronisation in the con- 
gested phase. The two trafSc patterns studied here are 
those of a single time deposition, and that of constant 
density traffic . Real life networks can have traffic pat- 
terns which wax and wane several times in a single day. 



However, synchronisation phenomena can be seen in real 
networks as well. We demonstrate this phenomenon in 
terms of the number of views at different web-sites, in 
the next section. 



VI. SYNCHRONIZATION FOR REAL TRAFFIC 
DATA 

We discuss phase synchronization in the Internet traf- 
fic data in the Indian Institute of Technology Madras, 
India. The data is collected for the number of views to 
different websites. The websites are www. google. com, 
www.gmail.com, www.yahoo.com, www.youtube.com 
and www.rediffmail.com. The data is counte d sp ecifi- 
cally for the given sites and not for subdomains [43 . Fig. 
I12f a) shows the total number of views for the five web- 
sites per day for a period of 92 days ,from 01/10/2008 to 
31/12/2008 {t in days on the x-axis). In Fig. \Mh) we 
plot the number of views per minute for 10*'' November 
2008 {t in minutes on the x-axis). 

As is evident from the plots the number of views for the 
website google is very large as compared to the rest. In 
Fig. [THit is seen that Nyiews for google show abrupt high 
peaks at ti = 13 (13*'' October), = 24 (24*'' October), 
h = 28 (28*'' October), i4 = 51 (20*'' November) and 
ts = 53 (22"'' November) (See gl]). 

In Fig. [T3]we plot the phase synchronization in num- 
ber of views for different websites. The phase is defined 
as in Eq.l. It is observed that the phase locking con- 
dition holds for the pairs (yahoo, gmail) and (youtube, 
rediffmail) (See Fig. fTSTa)). No such phase locking condi- 
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(c) t (d) t 



FIG. 11: (Color online) (a) Phase synchronization for constant density traffic for pairs of hubs as labeled and (b) Plot of order 
parameter r and average phase tp as a, function of time for the baseline mechanism. Messages are fed continuously at every 120 
time steps for a total run time of 12000. We consider 100 messages in a 100 x 100 lattice with 1% hub density and Z)st=142. 
(c) Plot of phase synchronization in queue length of the top most congested nodes ranked by their CBC values for (i) a = 0.05 
and /3 — 0.05; (ii) a = 0.4 and /3 = 0.4. = 500 points are chosen randomly in the lattice and connected by Waxman 
algorithm. We drop Nm = 100 messages on the lattice at 120 time steps, (d) The plot of global synchronization parameter r 
and ^ as a function of time t for a = 0.05 and /3 — 0.05 ; a — 0.4 and /3 = 0.4 for constant density traffic. 



tion exists between the pairs (google, yahoo) and (google, 
rediffmail) (See Fig. [TSTc)). This is due to two facts. 
First the number of views for www .google. com are much 
higher than those of the other websites. Secondly, the 
presence of abrupt peaks for google leads to larger fluctu- 
ations. The plot of global synchronization parameters r 
and ^' shows that the websites are synchronized in terms 
of number of views (See Fig. [TWb. d)). Larger fluctua- 
tions in r and are seen when all the five websites are 
taken into account (See Fig. [TSTb)). The fluctuations arc 
reduced when the website www .google.com is not taken 



into account (See Fig. fTSFd)). 

We observe similar behavior when the data is studied 
for 10*'* Nov 2008. It is observed that the two websites 
(youtube, rediffmail) peak together during the time in- 
terval t = 400 - 600 (Fig. [Hlja)). During this time 
interval the two websites are synchronized as shown in 
Fig. [TWa(ii)). As soon as the values of N^ews start de- 
creasing for both the sites, phase synchronization is lost. 
When we compared the sites google and yahoo it was 
observed that for yahoo, the number of views increases 
intermittently. Hence intermittent behavior of phase syn- 
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FIG. 12: (Color online) The plot of number views for a given website per day is obtained from proxy2.iitm.ac.in for a period 
of 92 days, from 01/10/2008 to 31/12/2008 (denoted by t in x-axis). The websites are www. google.com, www. gmail.com, 
www.yahoo.com, www .youtube.com and www. redif fmail.com. The data is counted specifically for the given sites and not for 
sub domains. 



chronization is observed for yahoo and google (Sec Fig. 
[TW b)). No phase synchronization is observed between 
google and redifFmail. This is similar to the absence of 
phase synchronization in queue lengths between higher 
CBC hubs and hubs of low CBC values as in Fig. d] for 
the 2 — d communication network. Also we observed that 
hubs of comparable CBC phase synchronize. Similarly 
in the Internet data, websites of comparable volume of 
traffic phase synchronize. 

We also study the traffic data of the global top five 
websites ranked by the percentage of global Internet 
users who visit the respective website. Fig. [Tbla] shows 
the plot of percentage of global Internet users, /„iet„s for 
the top ranked websites for a period of 16 days from 
08/02/2009 to 24/02/2009. The data has been collected 
from the website www. alexa.com. The websites arc 
www.google.com, www.yahoo.com, www.youtube.com, 
www.live.com and www.msn.com. Here also phase lock- 
ing is observed for the pairs of websites (Fig. [TSTb)). 



VII. CONCLUSIONS 

To summarize, we have established a connection be- 
tween synchronization and congestion in the case of two 
communication networks based on 2 — d geometries, a lo- 
cally clustered network, and a network based on random 
graphs, viz. the Waxman topology network. 

We first considered the case where many messages are 
deposited simultaneously on the lattice. We observed 
that the queue lengths of the top five hubs get phase 
synchronized when the system is in the congested phase 



i.e. the queue lengths at the hubs start piling up. Phase 
synchronization is lost when messages start getting de- 
livered to their destinations and queue lengths start de- 
creasing. Complete synchronization in queue lengths be- 
tween certain pairs of top five hubs is observed when 
decongestion strategies are implemented by connecting 
the top five hubs by gradient connections or by assorta- 
tive linkages. Phase synchronization is also seen between 
nodes of comparable CBC in the Waxman topology de- 
spite the fact that the degree distribution of the Waxman 
graphs is entirely different from that of the locally clus- 
tered network[50|. The phenomenon of synchronization 
in queue length in the congested phase is thus a robust 
one. We also observed that the phase synchronization 
does not exist between hubs of widely separated CBC-s, 
as can be seen from the lack of synchronization between 
the high CBC hubs and any randomly chosen hub (which 
turn out to have lower CBC values). 

It is has been seen, for typical configurations, that 
the central most region of the lattice is the most con- 
gested due to flow of messages from both sides of the 
latticefi^. The top hubs ranked by their CBC values 
were also seen to be located at the central region of the 
lattice and were close to each other. In such cases, the 
hub with the highest CBC value gets congested first and 
the remaining hubs follow, resulting in increase in queue 
lengths at the remaining hubs after a time lag result- 
ing in a phase difference between the queues at different 
hubs. This phase difference tends to a constant in the 
congested phase. Thus, phase synchronization in queue 
lengths exists pairwise for hubs of comparable CBC val- 
ues. Global synchronization is also seen between hubs of 
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FIG. 13: (Color online) (a) Phase synchronization in number of views for the pairs of websites (yahoo, gmail) and (youtube, 
rediffmail). (c) No phase synchronization is observed for (yahoo, google) and (google, rediffmail) as indicated by the large 
fluctuations in A^. (b) Large fluctuations in global synchronization parameter r and are seen when all the five websites are 
taken into account, (d) The fluctuations are reduced when the website www. google.com is not taken into account. 



comparable CBC-s. This is demonstrated by the behav- 
ior of the global synchronization parameters for the top 
five hubs. 

Similar phase synchronized behavior is observed when 
messages are fed on the lattice at a constant rate. The 
parameters of synchronization, such as the phase differ- 
ence and the global synchronization parameters, show 
small fluctuations in this case, indicating stronger syn- 
chronization as compared to simultaneous deposition. 

Thus, synchronization is associated with the inefficient 
phase of the system. Similar phenomena can be found 
in the context of neurophysiological systems (4l| and in 
computer networks [i^l ■ We observed that the most con- 
gested hub drives the rest. As soon as this hub decon- 
gests, the synchronization is lost. In the case of our com- 



munication network, this is usually the hub of highest 
CBC. Due to the master-slave relation between the most 
congested hub and the rest, there is cascading effect by 
which successive pairs lose synchronization. Cascading 
effects have been seen in other systems such as power 
grids [i^ and the Internet [i^]- It will be interesting 
to see if synchronization effects can be observed in these 
contexts. 

Finally we studied the Internet traffic data obtained in 
Indian Institute of Technology, Madras, India as well as 
the global Internet data obtained from the Alexa web- 
site. Despite the irregular pattern of traffic in this data, 
phase synchronization in the number of views is observed 
between two websites of comparable volume of traffic. 
Phase synchronization breaks down if the volume of traf- 
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FIG. 14: (Color online) (a(i)) The plot of number of views per minute (denoted by t in x-axis) for lO"^ Nov 2008. The websites 
are www .youtube.com and www.rediffmail.com. (a(ii)) Phase synchronization in number of views is observed for the pair 
(youtube, rediffmail) in the time interval t = 400 — 600. (b) The websites shown are www.google.com and www.yahoo.com. 
We observe intermittent phase synchronization in the number of views for (yahoo, google). The data is counted specifically for 
the given sites and not for sub domains. 
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FIG. 15: (Color online) (a) The percentage of global Internet users {f views) who visit the respective websites listed for a 
period of 16 days from 08/02/2009 to 24/02/2009 (denoted by t in x-axis). The websites are labeled as 1 {www.google.com), 
2 {www .yahoo.com) , 3 {www .youtube.com) , 4 {www .live.com) and 5 {www .msn.com) , where the numbers denote the global 
rank of the websites, (b) Phase synchronization in f views is observed for pairs of websites. 



fie changes abruptly. Global synehronization is also ob- 
served for the Internet data. This data also shows the 
existence of intermittent synchronization. 

These observations can be of great utility in the prac- 
tical situations. For example, e.g. synchronization can 
be considered as a predictor of congestion. Synehroniza- 
tion can also be used to detect changes in the pattern 
of traffic or to detect abnormal traffic from a given hub. 



The hub from where the attack originates can be easily 
identified via a synchronization effect. Synchronization 
in transport may also provide information about the way 
in which the network is connected. Thus our study may 
prove to be useful in a number of application contexts. 
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gested phase, whereas the S{t) — s are of order 1. Thus 
5{t) « q{t). 

[47] The value of r is seen to increase at the end of the decon- 
gestion phase. This is due to the fact that for the gradient 
mechanisms all queues are cleared and thus take the value 
zero at the end of the run. For the baseline mechanism 
the queues of Ath and 5th hubs are cleared while rest are 
trapped leading to a constant value of r which is less than 
one at the end of the run. 

[48] The data is obtained from the log files (generated by the 
SQUID software) of the proxy server in the IIT Madras 
campus 

[49] 13*'' October, 24*'' October, 20*'' November and 22'"* 
November were all dates for semester examinations in 
the Indian Institute of Technology Madras, India and 



students tend to access google more for tutorials and so- 
lutions available in the web, at these times. Again 28*'' 
October is the festival of Diwali which is a national hol- 
iday in India. The Internet users on campus appear to 
have spent most of this holiday browsing. The value of 
Nviews for all the websites reach their peak during the 
day time but decreases during night, contrary to the no- 
tion that web browsing increases during the night. This 
is due to the fact that Internet is unavailable between 
20:00 hours to 04:00 hours in the student hostels. 
[50] the geographically clustered network has a bimodal de- 
gree distribution [Ij], whereas the Waxman graph has de- 
gree distributions which depend on the values of a and 
/3H. 



